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Abstract. The group classification of variable coefficient quasilinear reaction- 
diffusion equations ut = Uxx + h{x)B{u) is carried out exliaustively. This 
became possible due to usage of a conditional equivalence group found in the 
course of the study of admissible point transformation within the class. 



1. Introduction 

Geometrical study of differential equations (DEs) has a long and distinguished 
history dating back to the second part of XIX century when the pioneering works 
of Gaston Darboux, Sophus Lie and Elie Cartan were published. Their ideas be- 
came a source for a number of developments including the theory of completely 
integrable systems and study of conservation laws. A brief but very nice review of 
both classical and modern treatments of geometrical study of differential equations 
is presented in [1]. One of the classical topics of such studies is Lie (point) sym- 
metries of differential equations. Lie proved that knowledge of continuous group of 
nondegenerate point transformations that leave an equation invariant allows one to 
reduce this equation to one with fewer independent variables. In many cases, this 
reduction procedure results in construction of a group-invariant solution in closed 
form. It is worthy to say that Lie symmetries represent by themselves powerful 
tool for finding exact solutions for partial differential equations (PDEs), and this 
is one of the most successful applications of geometrical studies of DEs [2, 3]. An- 
other feature of Lie symmetries is that they reveal equations which are important 
for applications among wide set of admissible ones. Indeed, all basic equations 
of mathematical physics, e.g., the equations of Newton, Laplace, Euler-Lagrange, 
d'Alembert, Lame, Hamilton-Jacobi, Maxwell, Schrodinger etc., have rich symme- 
try properties [4]. This property distinguishes these equations from other PDEs. 
Therefore, an important problem arises to single out from a given class of PDEs 
those admitting Lie symmetry algebra of the maximally possible dimension. This 
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problem is called group classification problem and is formulated as follows [3, 5]: 
given a class of PDEs, to classify all possible cases of extension of Lie invariance 
algebras of such equations with respect to the equivalence group of the class. 

At this stage the modern group analysis provides us with two main approaches 
for solving group classification problems. The first is algebraic one, based on sub- 
group analysis of the corresponding equivalence group [6, 7]. It results in complete 
group classification only if the class under study is normalized [8] . Roughly speak- 
ing, the class is normalized if any point transformation between two fixed equations 
from the class is induced by transformation from its equivalence group. The second 
approach is based on direct integration of determining equations implied by the 
infinitesimal invariance criterion. It is efficient usually only for classes of simple 
structure having one or two arbitrary elements of a single variable. Obviously, 
normalized classes or ones of simple structure do not exhaust a set of classes of 
PDEs or their systems which are important for application. To solve more group 
classification problems, a number of notions was introduced recently as well as new 
approaches were developed. These are, e.g., notions of admissible [8] (synonym: 
form-preserving [9]) transformations, generalized [10] and extended [11] equiva- 
lence groups, normalized class of DEs [8], equivalence groupoid [12], contractions 
of equations and conservation laws [13, 14]. Among new approaches it is wor- 
thy to mention the method of furcate split [15], the method of mapping between 
classes [16], the partition of a class into normalized subclasses [7, 8], etc. 

In this paper we solve the group classification problem for the class of variable 
coefficient semilinear reaction-diffusion equations of the form 



where h = h{x) and B = B{u) are arbitrary smooth functions of their variables, 
hBuu 7^ 0. Linear equations singled out of class (1.1) by the condition By^^ = 
are excluded from consideration since group classification of all second-order linear 
PDEs in two dimensions was performed by Lie (see [17]). 

Eqiiations from this class are used to model various phenomena such as mi- 
crowave heating, problems in population genetics, etc. (see, e.g., [18] and references 
therein). Theorems on existence and uniqueness of bounded solutions for equations 
of the more general form ut = u^x + F{t, x, u) were proven in [19]. 

Lie symmetries of certain subclasses of (1.1) are known. The group classifi- 
cation of constant coefficient equations from class (1.1) were studied by Dorod- 
nitsyn [20] (the results are adduced in handbook [5]). Class (1.1) includes the 
generalized Huxley equations 



whose Lie symmetries were studied in [21, 22]. There exists also certain intersection 
with the results on group classification of the classes 



(1.1) 



Ut = Uxx + h{x)B{u), 



Ut = Uxx + h{x)u'^{l - u), 



(1.2) 
(1.3) 



Ut = Uxx + H{x)u"' + F{x)u, H 0, 

ut = uxx + H{x)u^ +G{x), H^O, 
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which were obtained in [16]. Note that the group classifications for the general 
class of (l+l)-dimensional second-order quasilinear evolution equations 

ut = F{t, X, u, Ux)uxx + G{t, X, u, Ux), F ^ 0, 

was carried out in [6] . Nevertheless those results obtained up to a very wide equiv- 
alence group seem to be inconvenient to derive group classification for class (1.1). 

The structure of this paper is as follows. In the next section we investigate 
equivalence transformations in class (1.1). Lie symmetries are classified in Section 3. 

2. Equivalence transformations 

To solve a group classification problem for a class of differential equations it 
is important to describe point transformations that preserve the differential struc- 
ture of the class and transform only its arbitrary elements. Such transformations 
are called equivalence transformations and form a group [3]. By Ovsiannikov, the 
equivalence group consists of the nondcgcncratc point transformations of the inde- 
pendent and dependent variables and of the arbitrary elements of the class, where 
transformations for independent and dependent variables do not involve arbitrary 
elements, i.e., they are projectible on the space of independent and dependent vari- 
ables. If this restriction is neglected, then equivalence group is called a generalized 
equivalence group [10]. If new arbitrary elements appear to depend on old ones in 
some nonpoint (possibly, nonlocal) way, then the corresponding equivalence group 
is called extended [11]. Admissible transformation is a triple consisting of two fixed 
equations from a class and a point transformation that links these equations. The 
set of admissible transformations of a class naturally possesses the groupoid struc- 
ture with respect to the composition of transformations, and hence it is called the 
equivalence groupoid of the class [12]. We look for admissible transformations for 
class (1.1) using the direct method [9, 14]. 

Consider a pair of equations from the class under consideration, i.e., equa- 
tion (1.1) and the equation 

(2.1) ui = Uxx + h{x)B{u), 

and assume that they are connected via a point transformation T of the general 

form 

(2.2) i=T{t), i = X{t,x), u=U{t,x,u), 

where TtX^Uu ^ 0. We can restrict ourselves by the transformations of this form 
instead of general transformations t = T(t,x,u), x = X(t,x,u), u = U{t,x,u). 
It is due to the fact that class (1.1) is a subclass of more general class of (1-1-1)- 
dimensional quasi-linear evolution equations Ut = F{t, x, u)uxx + G(t, x, u, u^) with 
F 0, for which admissible transformations are proven to be of the form (2.2) [15]. 

We have to derive the determining equations for the functions T, X and U and 
then to solve them. Simultaneously we have to find a connection between arbitrary 
elements of equations (1.1) and (2.1). Substituting the expressions for the new 
(tildcd) variables into (2.1), we obtain an equation in the old (untilded) variables. 
It should be an identity on the manifold £ determined by (1.1) in the second-order 
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jet space with the independent variables {t, x) and the dependent variable u. 
To involve the constraint between variables of on the manifold C, we substitute 
the expression of tit implied by equation (1.1). The splitting of this identity with 
respect to the derivatives Uxx and Ux implies the determining equations for the 
functions T, X and U : 

(2.3) = 0, Xj" = Tu 2-^^" - ■ 



Uu Tt Xx 

(2.4) TthB - UuhB = Ut-Uxx-^Ux. 
Solving at first equations (2.3) we get that Tt > and 

ttI /■U\ 1 ^Ttt 2 ^ 

U = CU) exp —--—X — -—=x 
^ ^\ 8Tt 2 . 

where ?7°, a and C, are arbitrary smooth functions of their variables; £ = ±1. 
Then (2.4) can be written as 

(2.5) TthB - U'hB = J2(ui- - IjrUix - e^K^ «^ 

where = vP — 1. Investigating (2.5) for varying arbitrary elements h and B, 
we derive at first the usual equivalence group of class (1.1). 

Theorem 2.1. The usual equivalence group G~ of class (1.1) consists of the 
transformations 

A 

where 6j, j = 0, . . . ,5, are arbitrary constants with 6o5i64 ^ 0. 

It appears that there exist point transformations between equations from (1.1) 

which do not belong to C and form a conditional equivalence group. Moreover, 
this group is not usual but a generalized extended one. 

Theorem 2.2. The generalized extended equivalence group of the subclass 

(2.6) ut = u,, + /i(a;)(e"" + r) 
of class (1.1) is formed by the transformations 

t = Sft + 62, X = SiX + 63, U = S4U + Lp{x), 



t = 6it + 52, x = 6ix + 63, u = S4U + 65, h=-zj—h, B = SqB, 



h=-p;e ^i^h, h=—, f = e^i^\r- 

of 04 \ 04/1 

where r and 5j, j = 1, ... ,4, are arbitrary constants with SiS^ ^ 0. The transfor- 
mation component for r can be interpreted as the constraint for ip, 

fxx = hhir - fe~^'^). 
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Theorem 2.2 implies that class (2.6) reduces to the class 

(2.7) lit = uxx + h{x)e''^ 

by the transformation 

i=t, X = X, u = u + ip{x), 

where h{x) = e~'^'^^^h{x) and tpxx = rh{x). Class (2.6) is normalized. Therefore, 
the equivalence group of class (2.6) with r = can be found setting r = r = in 

transformations from the group G^p. 

Corollary 2.1. The usual equivalence group of the class 

ut = + h{x)e"'' 

consists of the transformations 

i=5ft + d2, X = Six + S3, u = S4U + S5X + Sq, 



Si _r^(^Sx+5e)u ~ 
n. = -^e " ^ "'ft, n = —. 

Sf ' d4 

where (5^ , = 1, . . . , 6, are arbitrary constants with SiS^ ^ 0. 



In the course of the study of Lie symmetries we will use the derived equivalence 
transformations for the simplification of calculations and for presenting the final 
results in a concise form. 

3. Lie symmetries 

We study Lie symmetries of equations from class (1.1) using the classical ap- 
proach [3] in combination with the method of furcate split [15]. We search for 
vector fields of the form 

Q = T{t, X, u)dt + ^{t, x, u)dx + r]{t, x, u)du 

that generate one-parameter Lie symmetry groups of a fixed equation C from 
class (1.1). These vector fields form the maximal Lie invariance algebra A™-^^ = 
^max^£-j q£ ^YiQ equation C. Any Lie symmetry generator Q satisfies the infinites- 
imal invariance criterion, i.e., the action of the second prolongation Q^^' of Q on 
the equation C results in the condition identically satisfied for all solutions of C. 
Namely, we require 



(3.1) Q^'''^{ut-Uxx-h{x)B{u)) 



= 0. 

c 



After elimination of Ut by means of (1.1), equation (3.1) can be regarded as a 
polynomial in the variables Ux, Uxx and Utx- The coefiicients of different powers of 
these variables should be zeros. This results in the determining equations for the 
coefficients r, ^ and rj. Solving these equations implies that r = r(f) and ^ = ^{t, x), 
which agrees with the general results on point transformations between evolution 
equations [9] . The remaining determining equations have the form 

(3.2) = Tt, r]uu = 0, 2rixu = £,xx - ^t, 

(3.3) r]hBu = {-^hx + (r?« - n) h)B + r]t- r]xx- 



6 



VANEEVA AND ZHALIJ 



Integrating equations (3.2) we get the following expressions for ^ and ry: 

£, = \TtX + a{t), rj = {-^Tttx^ ~ ^Utx + C{t)) u + T]°{t,x), 

where a, ^ and r]^ are arbitrary smooth functions of their variables. Then equa- 
tion (3.3) becomes 

(3 4) ((s'^"^^ + 5^*^ - C)w - r?°) hBu = {{Intx'^ + ^atx -C,+rt)h + 

+ {\TtX + (j) K) B + {Inttx^ + \crttx - Ct - jTtt) u-rj^ + rg^. 

It is called a classifying equation and should be solved simultaneously with respect 
to remaining uncertainties in the coefficients of infinitesimal generator Q, i.e., the 
functions r, cr, Q and rp , and arbitrary elements of the class, namely, h and B. 

In order to find the common part of Lie symmetries for all equations from 
class (1.1), we split with respect to the arbitrary elements in equation (3.4). This 
results in Tt = ^ = 7/ = 0. 

Proposition 3.1. The intersection of the maximal Lie invariance algebras of 
equations from class (1.1) (called kernel algebra) is the one- dimensional algebra 

A" = (dt). 

The next step is to classify possible extensions of A'^ using the method of furcate 
split [14, 15]. For any operator Q from A™*^^ the substitution of its coefficients 
into (3.4) gives some equation on B of the general form 

(3.5) {au + b)Bu = pB + qu + r, 

where a, b, p, q and r are constants which are defined up to nonzero multiplier. 
The set V of values of the coefficient tuple (a, b,p, q, r) obtained by varying of an 
operator from A™^^ is a linear space. Note that (a, 6) ^ (0, 0) and the dimension 
k = k{A"^'^^) of the space V is not greater than 2. Otherwise the corresponding 
equations imply that either B is linear in u or the system is incompatible. The 
value of k is an invariant of the transformations from C^. Therefore, there exist 
three G'"-inequivalent cases for the value of fc: A; = 0, A; = 1 and k = 2. We consider 
these possibilities separately. 

I. The condition fc = means that (3.4) is not an equation with respect to B but 
an identity. Therefore, B is not constrained and tu = <yt = C, = if = ^- We obtain 
that T = c\t + C2, a = Cs and the classifying equation on h 

(3.6) (icix + C3) + cih = 0, 

where Ci, C2, and C3 are arbitrary constants. If h is arbitrary we get the kernel 

algebra A'^ presented by Case of Table 1. It follows from (3.6) that the extensions 
are possible in two cases: either h = d{x+j3)^^ or h = 6, where jS and 6 are arbitrary 
constants, 6^0. Up to G~-equivalence 13 can be set to zero value and S to ±1 
depending on its sign. These two cases are presented by Cases 1 and 2 of Table 1. 

II. If = 1 then we have, up to a nonzero multiplier, exactly one equation of 
the form (3.5) with respect to the function B. Integration of this equation up to 
G"^ -equivalence gives three cases 

l.B = u"'+/3iu + p2, 2. B = e"" + Piu + 132, 3. B =ulnu + 02, 
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where /32 and m are arbitrary constants, m 7^ 0, 1. The next step is to substitute 
each of the three derived forms of B to (3.4) and subsequently to split it with respect 
to linearly independent functions of the variable u. Consider these cases separately. 

11. 1. Let B = M™ + (3iu + (32- If m 7^ 2 then the splitting equation (3.4) leads to 
the condition rf = Q and the three classifying equations on the function h 

(Inx + a)h^ + ((1 - n)W + Tt)h = 0, 02 {{^nx + a)h^ + {W + n) h) = 0, 
(1 - /3i) {^Ttx + a)h^ + {W+{l- pi)Tt) h-Wt + \Ttt = 0, 

where W = ^tux^ + \(JtX — Wt = The studying of compatibility of this 

system shows that if /32 7^ then t = cit + C2, a = c^, C = and the classifying 
equation on h is of the form (3.6). Therefore, this is subcase of the case k = Q 
considered above. 

If m 7^ 2 and ^2 = ^ then class (1.1) takes the form (1.2) with H{x) = h{x) and 
F{x) = j3\h{x). Using the results of [16] we derive that extension of is possible 
if /3i = and the function h take one of the forms presented by Cases 3, 4 and 5 of 
Table 1 . If /?i 7^ then the cases of extension of A'^ are subcases of Cases 1 and 2 
of Table 1. 

If B = V? + Piu + (32, then (3\ can always be set to zero by the translation 
u n- u — /3i/2 and then class (1.1) coincides with (1.3), where H{x?) = h{x) and 
G{x) = /32h{x). Using the results derived in [16] we obtain that the cases of Lie 
symmetry extension are either presented by Cases 1 and 2 of Table 1 if /32 7^ or 
by Cases 3-5 of Table 1 if = 0. 

11. 2. If i? = e" + /3iu + /32, then the classifying equations imply tu = <7t = C = 0. 
We obtain that t = cit + C2, a = C3 and the system of equations on h 

{{^cix + C3) hx + cih) /3i = 0, r/°/i + {^cix + C3) + cih = 0, 

((icia; + C3) /ix + cih) p2 + V°hpi + vL ' Vt = 0- 

If 7^ 0, then T]^ = and the classifying condition on h is exactly equation (3.6). 
This case can be included to the general case with arbitrary B. If (3^ = 0, then 
/32 can be set to zero using the transformation from the conditional equivalence 
group G^p. Then ry" = c^x + C5, where C4 and C5 are arbitrary constants and the 
remaining classifying equation on h takes the form 

(3.7) (^cia; + C3) + {ciX + C5 + ci)h = 0. 

Combined with the multiplication by a nonzero constant, each transformation from 
the equivalence group is extended to the coefficient tuple of the above equation 
in the following way 

Ci=KCi, C3 = K (C3(5i - ici^s) , C4 = ^ (C4 + ici^s) , 

Ol 

h = Y (^5^1 + C3(5i(55 - C4^3 - |ci (53(55) . 

Here k is an arbitrary nonzero constant. Using these transformations we derive the 
following statement. 
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Table 1. The group classification of tlie class ut — Uxx + h{x)B(u), hB^u 7^ 0. 



no. 


B(u) 




Basis of A™"" 





V 


V 


a* 


1 


V 




a, , 2/:a, + xa. 


2 






a, . a,. 


3 




Sx" 


dt, 2(m - l)tat + (to - l)xd^ - (« + 2)u9„ 


4 


u" 




St, (1 - ■m)dx + udu 


5 


u" 


<5 


dt, dx, 2 (to - l)tdt + (to - - 2u9„ 


6 


e" 




at, 2ta( + - (s + 2)a„ 


7 


e" 




at, a^ T 2xdu 


8 


e" 




dt, dx, ^tdt + xdx - 2du 


9 


ulnu 




dt,dx, e^^udu, e^\dx - ^xudu) 



Here 5, m and s are arbitrary constants, m ^ 0, 1, s ^ 0, 5 = ±1 mod C 

Lemma 3.1. Up to C^^^- equivalence the parameter tuple (01,03,04,05) can be 
assumed to belong to the set {(1, 0, 0, 05), (0, 1, ±2, 0), (0, 1, 0, 0), }, where 05 is an 
arbitrary constant. 

Integration of (3.7) up to G^p-equivalence leads to the forms of h presented 

in Cases 6 8 of Table 1. In Case 6 wc use the notation s = —2(65 + 1). 

II. 3. If _B = ulnu + l32, then 77*^ = and the remaining classifying conditions are 

(irtx + a) ha; + Tth = 0, (^Tttx^ + ^atx - () h = \TtttX^ + \<7ttx - Ct - jnt, 

p2 {{\Ttx + (j)hx + (gTttX^ + \atx -C + Tt)h)= 0. 

Investigation of this system implies that if /32 7^ or /32 = but hx ^ Q, then 
Ttt = cTt = C = 0- Therefore t = cit + 02, tr = 03 and the classifying equation on 
h is (3.6). I.e. we get nothing but subcases of Cases 1 and 2 of Table 1. But if 
/32 = and h = 6 = const, then the classifying equations leads to the conditions 
Tf = att-Sat = Ct-SC = 0. Therefore r = ci, ^ = 02+030**, 77 = (-fo3X + 04) e^*u 
and we have extension of A'~' on three Lie symmetry operators (Case 9 of Table 1). 



III. Let k = 2. We choose a basis {{ai,bi,Pi,qi,ri), i = 1,2} of the space V of 
tuples {a,b,p,q,r) associated with A™'^^. The system on B is of the form 

(oiu + bi)Bu = PiB + qiu + n, 

{a2U + b2)Bu = P2B + q2U + r2- 

The determinant of the matrix M = ( I is nonzero, since otherwise B is 

\a2 b2 J 

linear in u. Therefore the system can be rewritten in the form 

uBu =p'iB + q[u + r[ , 
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Bu =P2B + q'2U + r'2, 

- fe) • (:0 - (:;) • (.1) - ft) ■ - 

the inverse matrix of the matrix M. Studying the latter system we get that it is 
compatible only if B is either linear or quadratic. The case B is quadratic in u 
is considered already in the course of the study of the case B = + ^lU + ^2- 
So, the case k = 2 docs not lead to any new case of Lie symmetry extension for 
equations (1.1) with Buu 7^ 0. 

The group classification problem for class (1.1) is solved exhaustively. The 
results are summarized in Tabic 1. It is important to note that group classification 
of subclass (2.6) is carried out up to the G^p-equivalence, whereas all other cases 
are classified up to the usual G'"-equivalence. 

4. Conclusion 

In this paper we solve the group classification problem for the class of (1+1)- 
dimensional quasilinear diffusion equations with a variable coefficient nonlinear 
source (1.1) which arise as mathematical models in problems of mathematical biol- 
ogy and other applied areas [18]. The results of group classification can be applied 
for searching closed form solutions via the classical reduction method. The knowl- 
edge of Lie symmetries is also necessary for finding nonclassical symmetries (called 
also Q- conditional symmetries or reduction operators) of equations (1.1). This will 
be a subject of a forthcoming paper. 

Acknowledgement. The authors are grateful to Roman Popovych and Iryna 
Yehorchenko for useful discussions and constructive suggestions. OV thanks the 
Organizing Committee of the XVII Geometrical Seminar for the hospitality and 
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